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Onsager reciprocity in premelting solids 
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The diffusive motion of colloidal particles dispersed in a premelting solid is analyzed within the 
framework of irreversible thermodynamics. We determine the mass diffusion coefficient, thermal 
diffusion coefficient and Soret coefficient of the particles in the dilute limit, and find good agreement 
with experimental data. In contrast to liquid suspensions, the unique nature of premelting solids 
allows us to derive an expression for the Dufour coefScient and independently verify the Onsager 
reciprocal relation coupling diffusion to the flow of heat. 
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The melting of any material is normally initiated at 
one of its free surfaces at temperatures below the bulk 
melting temperature, T^, by the formation of a thin liq- 
uid - interfacially premelted - film. This surface phase 
transition has been observed at the interfaces of solid 
rare gases, quantum solids, metals, semiconductors and 
molecular solids including ice, allowing the liquid phase 
to persist in the solid region of the bulk phase diagram 
[1]. When T « T„i the premelted film is thicker than the 
correlation lengths of the solid- liquid interfaces and hence 
the total free energy of the (planar-planar) system is rep- 
resented as a linear combination of bulk and interfacial 
terms (this latter contains both the interfacial energy and 
the long ranged volume-volume interactions which are 
also forces/area in such a system). The total free energy 
of a curved interfacially melted system includes in addi- 
tion the Gibbs-Thomson effect, which, while contributing 
to the premelted film thickness d, one can prove produces 
no net thermodynamic force over a closed surface [2!|. 
Nonetheless, when a premelted film forms around a for- 
eign particle within a subfreezing solid, the particle can 
migrate under the influence of a temperature gradient, 
which produces a thermomolecular pressure gradient, a 
phenomenon referred to as regelation or thermodynamic 
buoyancy [2] [3] . Here we analyze this and related phe- 
nomena within the framework of irreversible thermody- 
namics and demonstrate that in addition to motion under 
the influence of a temperature gradient, the particle can 
undergo constrained Brownian motion owing to thermal 
fluctuations in the premelted film. We determine the 
diffusivity in the dilute limit and relate it to the Stokcs- 
Einstein diffusivity of particles in the bulk melt. Fur- 
thermore, the motion of a particle is shown to induce 
a reciprocal effect - a heat flux due to the release and 
absorption of latent heat on opposite sides of the parti- 
cle. We show that this effect is described by the Onsager 
relation coupling mass diffusion to heat flux. 

We consider a suspension of spherical particles of ra- 
dius R randomly distributed in a bulk solid. For tempera- 
tures T near but less than Tm , each particle is surrounded 
by a premelted film of thickness d that facilitates both 
the particles' constrained Brownian diffusion and their 
directed motion parallel to temperature gradients (figure 
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FIG. 1: Schematic diagram of a premelting solid containing 
colloidal inclusions at local volume fraction (h. 



[T|). From linear irreversible thermodynamics |3] the 
equations describing the fiux of particles and heat are 
(see Appendix) 

3 ^ ~DV4>+{l~(j))(j)DTyT, (I) 

(l=-kWT+{l~(l>)n^T^DgWcb, (2) 

where J = (j)V is the particle volume fiux, with (j) the 
volume fraction of particles and V the local average par- 
ticle velocity. In equation ([2| q is the heat fiux, k is 
the thermal conductivity and H^ is the derivative of the 
osmotic pressure of the particles with respect to volume 
fraction. The coefficients D, Dt and Dq are, respec- 
tively, the mass diffusion coefficient, thermal diffusion 
coefficient and Dufour coefficient. Here, in the dilute 
limit, we determine these quantities as functions of the 
undercooling AT = Tm — T . 

The second term on the right hand side of equation 
([T| accounts for the Soret effect, in which a temperature 
gradient generates motion of the particles. We define Dt 
to be positive when the particles migrate towards their 
warm sides. Whereas in liquid suspensions the Soret ef- 
fect is a complex phenomenon, with the Soret coefficient 
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St = Dt/D often changing sign as the temperature in- 
creases O [7], in premelting sohds the particles always 
appear to migrate towards higher temperatures. Further- 
more, the Soret coefficient is several orders of magnitude 
larger than in liquid systems. 

The second term in equation ([2| represents the Du- 
four effect, in which a concentration gradient gives rise 
to a flow of heat. Using the principle of microscopic re- 
versibility, Onsager proved that the Dufour coefficient 
Dq is equivalent to the thermal diffusion coefficient Dt 
[U[S]. There has, however, been significant controversy 
over the range of validity of the Onsager relations [HI IS] 
and ongoing attempts to prove them in particular cases 
[ini [m [12] . We show that, in the dilute limit, so long as 
the temperature gradient satisfies linear irreversible ther- 
modynamics and is constant on the scale of the particle, 
the Onsager relation is satisfied. 

Colloidal particles dispersed throughout a solid near 
its melting temperature undergo random motion owing 
to thermal fluctuations in their premelted films. Their 
random walk is driven by the thermodynamic force 



Fb = --vn, 



(3) 



where 11 is the osmotic pressure of the particles and Vp — 
|7ri?^ is the particle volume [T^]. For temperatures near 
Tm in the dilute limit the osmotic pressure is given by 
the van't Hoff equation 



If — kbTjm 



(4) 



where ki, is Boltzmann's constant. 

Owing to viscous flow in the premelted films, the par- 
ticle velocity V is accompanied by a net lubrication force 
given by 



F = —V 



(5) 



where K= d^/6R is a permeability coefficient [TT and rj 
is the dynamic viscosity of the premelted liquid film. 

In mechanical and thermal equilibrium, the total force 
F^ + Fb= 0, giving 

Comparing this result with equation ([T]) yields 

Hence we see that the premelting-controUed Brownian 
diffusivity D differs from the Stokes-Einstein diffusiv- 
ity Do of particles in bulk liquid by the factor 3d'^/4i?'^. 
Preliminary studies using X-ray photon correlation spec- 
troscopy show promise in testing equation ([t]) [15] . 

In the presence of a temperature gradient variations 
in the thickness of the premelted film over the surface of 



an isolated particle lead to a thermomolecular pressure 
gradient force causing the particle to move [ll[2j[3]. By 
integrating the fluid pressure over the surface of a particle 
this force has been obtained as 



Ft = Up'PtVT, 



(8) 



where Vt = Pslm/Tm is the thermomolecular pressure 
coefficient characterizing the magnitude of the inter- 
molecular forces responsible for the premelted film [l], 
and ps and Qm are the mass density and latent heat of 
fusion, respectively, of the bulk solid T . (We express Eq. 
([s]) differently than in for consistency with the present 
extension of that work.) In mechanical equilibrium the 
thermomolecular force is balanced by the viscous force 
F,; leading to the relation 



V=-7'tVT, (V0 = O). 
V 



(9) 



Comparing ([T|) with ^ in the limit 0-^0 yields an 
expression for the thermal diffusion coefficient 
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6Rt]T„, 



(10) 



This relation can be tested experimentally using the re- 
sults of a study in which a 12.7 /xm glass bead encased 
in ice was moved by imposition of a linear temperature 
gradient 16J. Figure[2]shows experimental measurements 
(stars) of Dt = V/VT as a function of the undercooling. 
Using measurements of d = d{AT) [TTJ [TB] in equation 
( 10 1 we can predict Dt{AT) which is shown by the solid 
lines. Clearly the agreement between theory and exper- 
iment is good. The deviation at small AT may be due 
to an overestimation of the film thickness by Ishizaki et 
al. |17j owing to curvature effects in their system, which 
become important as AT approaches 0. 
The Soret coefficient St is defined as 
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(11) 



Combining (jlj, Q and (10 1 with (111 gives 

St=^='-^, (12) 



Dt _ Vt 



in terms of fundamental quantities, where here we em- 
phasize that £t = kjyTm/vp is the ensemble averaged ki- 
netic energy/volume (thermal kinetic energy) of a Brow- 
nian particle. Thus in premelting solids the Soret coeffi- 
cient represents the ratio of the thermomolecular pressure 
coefficient to the thermal kinetic energy of a Brownian 
particle and is predicted to be independent of premelted 
film thickness and undercooling. For micron sized par- 
ticles in ice equation (12) yields St ~ Kf" thereby 



illustrating the dominance of the thermomolecular force 
over the Brownian force in premelting solids. As a com- 
parison, for /im sized particles in liquids St ^ 10~^ 
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FIG. 2: Thermal diffusion coefficient of silica particles in 
ice. The data (stars) are from Romkens and Miller [16) . 
while the theoretical curves are from equation ( 10 1 using d 



obtained from the averaged experimental measurements of 
Ishizaki et al. [17] and Engemann et al. [18] at each AT. 
The properties of ice and water at Tm are: ps = 920kg/m^, 
7? = 0.0017Ns/m^ = 3.34 x 10^ J/kg. 



expression for the concentration dependence of the Soret 
coefficient: 



(1 



(15) 



The Onsager relation will break down when the temper- 
ature gradient varies significantly over the surface of a 
particle because, as noted by Rempel et al. [2j, under 
these circumstances the thermomolecular pressure gradi- 
ent will vary over the surface of the particle and hence 
the expression ([8| for the thermomolecular force will no 
longer hold. Such a constraint on VT is independent of 
the issue of the validity of linear irreversible thermody- 
namics as embodied in equations ([T]) and This is 
because while VT may be sufficiently small for these re- 
lations to be satisfied it may still vary on the particle 
scale, say, for large particles. Nonetheless, for most ex- 
perimental situations envisaged the temperature gradient 
is sufficiently small that the Onsager relation holds. 

We have shown the complimentary roles that interfa- 
cial premelting plays in both the regelation and Brownian 
motion of colloidal particles in bulk solids. In the dilute 
limit, within the framework of linear irreversible thermo- 
dynamics, the Brownian difFusivity and Soret coefficient 
have been determined. Furthermore, in premelting solids 
the heat of transport can be explicitly calculated, lead- 
ing to a verification of the Onsager relation coupling mass 
diffusion to the flow of heat. 



[B]. It is hoped that these forms of St provide an ideal 
experimental target. One promising method for obtain- 
ing both D and simultaneously is thermal diffusion 
forced Rayleigh scattering [TO] . 

When foreign particles of volume Vp are inserted into a 
premelting host material the heat of fusion absorbed by 
the solid per unit volume of particles is ps^m- Whereas, 
when particles move through an initially isothermal solid 
they melt material ahead, absorbing latent heat, while 
freezing solid behind, releasing heat. The volumetric heat 
of transport characterizing this phenomenon is defined by 



VT=0 



= -(1 - <t>)Ii^Tr^ 



' D 



(13) 



where the second equality follows from ^ and (U. 
Hence, physically, Q* represents the sensible heat re- 
leased by an arbitrarily isolated part of the system per 
unit volume of particles diffusing into it (e.g., [201 f2T| ) 
and we obtain Q* = —psqm- 

We verify the reciprocal relation coupling mass diffu- 
sion to heat flux in our system by inserting (|4]) , ([t]) , ([To]) 
and Q*^ = —psq-m into equation ( 13 1 and taking the limit 
6 ^ Qto obtain 



(14) 



If the Onsager relation (14) holds also at higher concen- 
trations, (111, (13 1 and (14 1 can be used to obtain an 
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APPENDIX: MASS AND HEAT FLUXES 

Here we show that equations ([T| and (|2| can be de- 
rived from the irreversible thermodynamic treatment of 
the flux of mass and heat in a two-component hydrostatic 
mixture (cf., equations (XI.226) and (XI.227) of [4 ), viz.. 



Ji = —pDVci — pciC2Dt^T, 



3'^ = -kVT-p, 



dp-i 
dci 



TDqVci. 



(A.l) 



(A.2) 



T,P 



Here Ji — pi{^i — v) is the mass flux of component 
1 relative to the barycentric velocity v = civi + C2V2, 
Cj — Pj / p is the mass fraction of j which has a partial 
mass density pj, p — pi + p2 is the mixture density, and 
Vj is the velocity of j with respect to the laboratory 
frame. The so-called reduced heat flux is J^, where T is 
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absolute temperature, P is the mixture pressure and iij 
is the chemical potential of j. By definition 



Cl 



C2 = 1 



and 



Ji + J2 = 0. 



(A.3) 



With (A.3 1 equations (A.l) and (A. 2 1 can be written as 
J2 = -pDWc2 + pciC2DtVT, (A.4) 



dm 

dc2 



TD„Vc2. 



(A.5) 



T,P 



It is convenient to write equations (A.4 1 and ( A.5| in 
terms of the partial density p2- To achieve this we use 
the thermodynamic identity (equation XI. 99 of |4]) 



dp2 



= P vi, 



(A.6) 



where vi is the partial specific volume of component 
approximately equal to 



where ps is the density of pure component 1 at T and P. 
Additionally, we make use of the fact that the mass flux 
can be written in terms of the volume average velocity 
v'' = pivi^i + P2V2^2 using the thermodynamic identity 
p\V\ + /92'^2 = 1 to obtain 



J2 = pv\'i2 



(A., 



where J2 = P2(v2 — v"). 



this and (A.5) become 



With (A. 61, we rewrite (A.4| in the form of (A.8I and 



-D\I P2 + VxP\P2P>T^T, 



-fcVT- 



dp2 



TDqVp2 



(A.9) 



(A.IO) 



T,P 



For a dispersion of hard spherical particles 1^2 = 
where pp is the constant mass density of an individual 
particle, in which case (A.9) and (A.IO I can be written 
as 



-kVT 



Pi 



dpi 



(A.ll) 



(A.12) 



T,P 



(A 7) where J — 4>{v2 — v°) is the particle volume flux, — 
P2/ Pp is the volume fraction of particles and q = J^. As 



shown by de Groot and Mazur, in closed systems in which 
ui and V2 are constant, the volume average velocity Vg is 
zero, in which case ( A.lip reduces to ([ij. Inserting the 
thermodynamic relation (9/ii/9(/))t„,p„ — ^(1/Ps)ll0 
into ( |A.12 1 gives equation ([2]), where 1 — (p — p\l Ps and 
= (MI/90)t„,p„ at a reference state of bulk coexis- 
tence. 
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